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Abstract 

Background: Genonnic innprinting is an innportant epigenetic process involved in regulating placental and foetal 
growth. Imprinted genes are typically associated with differentially methylated regions (OMRs) whereby one of the 
two alleles is DNA methylated depending on the parent of origin. Identifying imprinted DMRs in humans is 
complicated by species- and tissue-specific differences in imprinting status and the presence of multiple regulatory 
regions associated with a particular gene, only some of which may be imprinted. In this study, we have taken 
advantage of the unbalanced parental genomic constitutions in triploidies to further characterize human DMRs 
associated with known imprinted genes and identify novel imprinted DMRs. 

Results: By comparing the promoter methylation status of over 14,000 genes in human placentas from ten 
diandries (extra paternal haploid set) and ten digynies (extra maternal haploid set) and using 6 complete 
hydatidiform moles (paternal origin) and ten chromosomally normal placentas for comparison, we identified 62 
genes with apparently imprinted DMRs (false discovery rate <0.1%). Of these 62 genes, 11 have been reported 
previously as DMRs that act as imprinting control regions, and the observed parental methylation patterns were 
concordant with those previously reported. We demonstrated that novel imprinted genes, such as FAM50B, as well 
as novel imprinted DMRs associated with known imprinted genes (for example, CDKNIC and RASGRFl) can be 
identified by using this approach. Furthermore, we have demonstrated how comparison of DNA methylation for 
known imprinted genes (for example, GNAS and CDKNIQ between placentas of different gestations and other 
somatic tissues (brain, kidney, muscle and blood) provides a detailed analysis of specific CpG sites associated with 
tissue-specific imprinting and gestational age-specific methylation. 

Conclusions: DNA methylation profiling of triploidies in different tissues and developmental ages can be a 
powerful and effective way to map and characterize imprinted regions in the genome. 



Background 

Genomic imprinting is a phenomenon in which one of 
the two alleles of a gene is expressed in a parent-of-ori- 
gin manner [1]. Imprinted genes are thought to be parti- 
cularly important to placental and foetal growth and 
development and may help regulate growth in response 
to maternal and foetal signals in utero [2]. To date, 
around 60 imprinted genes have been identified in 
humans (http://www.geneimprint.com), largely after first 
being identified in mice or through characterization of 
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specific imprinting disorders such as Prader-Willi syn- 
drome and Angelman syndrome or Beckwith-Wiede- 
mann syndrome. However, many genes are imprinted in 
mice but are not known to be in humans, for example. 
Impact [3]. Furthermore, many genes are imprinted only 
in specific tissues, for example, UbeSa, which is mater- 
nally expressed in the brain but biparentally expressed 
in other tissues [4], or may be polymorphically 
imprinted, for example, IGF2R [5]. These issues compli- 
cate the discovery and characterization of imprinted 
genes in humans. 

The importance of imprinted genes for placental and 
foetal development was initially demonstrated in mice 
by observations that parthenogenetic embryos (maternal 
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origin, digynic diploid) could show embryonic differen- 
tiation but failed to form extraembryonic components 
[6]. In contrast, androgenetic embryos (paternal origin, 
diandric diploid) had poorly developed embryos, but the 
trophoblasts showed extensive proliferation [7]. The par- 
allel observations in humans are ovarian teratomas 
(parthenogenetic), which are a rare form of tumour that 
consists of a variety of embryonic tissues or organs but 
no placental tissues, and complete hydatidiform moles 
(CHMs) (androgenetic), which consist of abnormal pla- 
cental growth characterized by trophoblast hyperplasia 
but no (or rare) embryonic structures. The parental con- 
flict theory developed to explain the evolution of 
imprinted genes [8] suggests that paternally expressed 
genes tend to promote growth of the offspring at the 
expense of the mother, while maternally expressed 
genes act as growth-limiting factors to conserve mater- 
nal resources [8]. 

Most imprinted genes possess differentially methylated 
regions (DMRs) whereby allelic methylation depends on 
the parent of origin [1]. DMRs established through the 
germline are called 'gametic' or 'primary' DMRs. These 
often coincide with imprinting control regions (ICRs), 
which regulate gene expression and further epigenetic 
modifications [9-11]. Their methylation status is thought 
to be maintained in all somatic lineages once acquired. 
Other DMRs, called 'somatic' or 'secondary' DMRs, are 
established after fertilization and may be tissue-specific 
[10,11]. 

Since most imprinted genes contain DMRs, comparing 
DNA methylation profiles between tissues with unba- 
lanced parental constitutions provides an approach to 
identify and characterize imprinted genes in the gen- 
ome. One approach is to compare the DNA methylation 
profile of paternally derived CHMs to that of maternally 
derived ovarian teratomas [12]. Indeed, several novel 
imprinted genes have been identified previously by using 
this strategy [13,14]. However, such comparisons are 
limited by the fact that the tissues present in ovarian 
teratomas and CHMs are highly abnormal and are not 
of comparable origin, with teratomas being embryonic 
and CHMs being strictly placental. Many differences 
may reflect tissue-specific methylated genes, since tis- 
sue-specific DMRs are numerous and are established in 
early pregnancy [15]. CHMs also present with highly 
proliferative trophoblasts that can lead to increased risk 
of choriocarcinoma, and hypermethylation of nonim- 
printed genes has been reported in CHMs [16]. 

In humans, triploidy (the presence of three complete 
haploid genomes) occurs spontaneously in 2% to 3% of 
pregnancies, and, while such pregnancies frequently end 
in miscarriage, they can survive into the foetal period 
and, very rarely, to term [17]. We propose that a com- 
parison between diandric and digynic triploidies, in 



which development is much less severely altered than in 
CHMs and teratomas, provides a powerful approach for 
the identification and characterization of imprinted 
genes in the human genome. The diandric triploid phe- 
notype (two paternal plus one maternal haploid gen- 
omes) is characterized by a normal-sized or only 
moderately growth-restricted foetus with a large and 
cystic placenta with trophoblast hyperplasia, while the 
digynic triploid phenotype (two maternal plus one pater- 
nal haploid genomes) is characterized by an intrauterine 
growth-restricted foetus and a very small placenta with 
no trophoblast hyperplasia [17]. Importantly, embryo 
and foetal development are largely similar between dia- 
ndric and digynic triploidy, with growth differences 
likely arising largely as a consequence of differences in 
placental function [18]. Furthermore, while small, digy- 
nic placentas have a grossly normal structure. Diandric 
placentas show features similar to a CHM, but their 
development is much less severely altered than in a 
CHM, and the placenta can support growth of a foetus 
at least to some degree. 

Although it was previously suggested that DNA 
methylation may be less important in regulating 
imprinting in placental tissue as compared to foetal tis- 
sue, we recently demonstrated that the DNA methyla- 
tion status of many known imprinted DMRs is strictly 
maintained in triploid placentas and can be used to dis- 
tinguish diandric from digynic triploidy [19]. Therefore, 
in the present study, we compared the DNA methyla- 
tion profiles of placentas from diandric and digynic tri- 
ploidies using the Infinium HumanMethylation27 
BeadChip array (lUumina, Inc., San Diego, CA, USA), 
which targets over 27,000 CpG loci within the proximal 
promoter regions of approximately 14,000 genes [20]. 
Methylation levels in chromosomally normal placentas, 
CHMs and maternal blood samples were used as refer- 
ence points for comparison. Using this strategy, we 
identified the majority of known imprinted ICRs on the 
array and many novel imprinted DMRs in the genome. 
For a subset of genes, we identified expressed poly- 
morphisms and informative mother-placenta pairs, 
which were used to demonstrate parent-of-origin biases 
in allelic expression. We also demonstrated that com- 
plex DNA methylation domains that regulate imprinted 
genes can be mapped by comparing the methylation 
patterns in different tissues and different gestational 
ages of placentas. 

Results 

DNA methylation profile analysis in placenta and blood 
samples 

To generate DNA methylation profiles from triploidies, we 
assayed placental DNA from ten diandric and ten digynic 
triploidies on the Illumina Infinium HumanMethylation27 
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BeadChip panel In addition, ten chromosomally normal 
placentas, 6 CHMs (diandric diploid, no maternal contri- 
bution) and ten maternal whole-blood samples were 
included for comparison. After background adjustment 
and normalization, we performed unsupervised hierarchi- 
cal clustering with all the samples based on a distance 
measure of 1-r, where r is the Pearson correlation coeffi- 
cient between different samples. This revealed three dis- 
tinct groups of clusters: CHMs, triploid and normal 
placentas, and blood (Figure 1). The blood cluster is more 
distant from the two other clusters of placentas, confirm- 
ing that there are many DNA methylation differences 
between blood and placenta [21-23]. Although CHMs are 
trophoblast- derived, they show a distinct methylation pro- 
file from the triploid and normal placentas, which prob- 
ably reflects not only the lack of a maternal genome but 
also the abnormal development of such tissue. Within tri- 
ploid and normal placentas, digynic and diandric triploid 
placentas are clearly separated by their methylation pro- 
files, but, interestingly, they are not separated from the 
chromosomally normal placentas (Figure 1). This suggests 
that methylation profiles of triploid placentas closely 
resemble those of chromosomally normal placentas, but 



that digynic and diandric triploid placentas have distin- 
guishing DNA methylation differences. 

Although clustering can be biased by gender differ- 
ences resulting from inactivation of an X chromosome 
in females (that is, higher methylation of the X chromo- 
some CpG islands in female than in male samples) 
[22,24], there is no preferential clustering of samples by 
gender within the triploid and normal placenta cluster 
(Figure SI A in Additional file 1). There is a small differ- 
ence in gestational age (about three weeks apart on 
average) between diandric and digynic placentas {P < 
0.01) (Table SI in Additional file 2), but this also cannot 
explain the distinct clustering patterns, since the gesta- 
tional ages of the two groups are largely overlapping 
(Table SI in Additional file 2). 

We further compared the average DNA methylation of 
probes between the five sample groups (digynic triploid 
placentas, diandric triploid placentas, normal placentas, 
CHMs and blood) (Figure SIB in Additional file 1). As 
expected, the correlation of average probe methylation 
values between different sample groups is consistent 
with that observed in the cluster analysis. In general, 
blood has the most distinct DNA methylation profile 
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Figure 1 Unsupervised clustering of triploid and normal placentas with CHMs and blood samples demonstrates that each tissue type 
has a distinct methylation profile. Sample names are shown with labelling of corresponding tissue types. Samples were clustered by 
hierarchical clustering of (3 values based on 1-r (lllumina GenomeStudio software), where r represents the correlation coefficient between 
samples. Digynic triploids are indicated by red boxes, diandric triploids by blue boxes and normal placentas by green boxes. 
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with a greater number of highly methylated probes (Fig- 
ure SIB in Additional file 1). Triploid and normal pla- 
centas are highly correlated with regard to their 
methylation profiles (r = 0.99), while CHMs are more 
similar to diandric and normal placentas (r = 0.98) than 
to digynic placentas (r = 0.96). 

Comparison of DNA methylation profiles between 
placentas from diandric and digynic triploidies 

After comparing methylation between diandric and digy- 
nic placentas by performing Student's ^-test for all 
probes, nearly 2,500 probes were identified with a P 
value < 0.01, which is nearly ten times more than 
expected by chance (Figure SIC in Additional file 1). To 
adjust for multiple testing and identify candidates with a 
very high likelihood of representing true differences, we 
used a stringent cutoff of <0.1% false discovery rate 
(FDR) by using the Significance Analysis of Microarrays 
(SAM) program with 1,000 permutation comparisons 
for each sample [25]. To further focus on the most 
meaningful differences, we also considered only probes 
with more than 15% absolute magnitude difference 
between the mean methylation of diandric and digynic 
triploidies. While we expected a theoretical difference of 
33.3% for imprinted sites, we used a lower cutoff 
because we have observed that the actual methylation 
difference may vary for some known imprinted genes 
[19] and that there may be biases in the lUumina array 
that result in a nonlinear relationship between the esti- 
mated methylation p value and actual methylation. In 
total, 122 probes were identified with <0.1% FDR and 
average absolute methylation difference >15% (average 
absolute A (3 >0.15 from the lUumina array). Probes 
with higher average methylation in diandric than digynic 
triploidies were designated putative paternal differen- 
tially methylated loci (DML), and probes with higher 
average methylation in digynic than diandric triploidies 
were designated putative maternal DML. Plotting DNA 
methylation of putative DML in all samples from dia- 
ndric against digynic triploidies showed a clear separa- 
tion of methylation values of paternal and maternal 
DML (Figure SID in Additional file 1), suggesting that 
most of the identified differentially methylated probes 
are consistently methylated within each sample group 
without much overlap as expected on the basis of our 
application of stringent statistical criteria. 

As some methylation differences between diandric and 
digynic triploids could theoretically arise as a result of 
secondary effects, such as altered cell composition, the 
validity of the identified putative imprinted DML was 
further evaluated by verifying that the methylation levels 
of diandric CHMs and chromosomally normal placentas 
fit the expected pattern (Figure 2). The average methyla- 
tion in CHMs was closer in value to that of diandric 



triploidies (Figures 2A and 2C), while that for normal 
placentas fell between that for diandric and digynic tri- 
ploidies for the majority of putative DML (Figures 2B 
and 2D) as would be expected for imprinted DMRs. The 
putative maternal DML were more strongly correlated 
with normal placentas than paternal DML, while puta- 
tive paternal DML tended to have higher correlation 
with CHMs than maternal DML (Figure 2 and Figure 
SIE in Additional file 1). CHMs showed particularly low 
correlation for maternal DML compared with other pla- 
cental groups, which was largely due to the low average 
methylation of putative maternal DML in CHMs as well 
as more variability in values for CHMs (Figure 2D). 

Fourteen probes failed to follow the expected relative 
methylation patterns between the between groups (nor- 
mal placentas with an average methylation level between 
that for diandric and digynic placentas and CMHs with 
an average methylation level closer to that in diandric 
placenta), and these loci were eliminated as candidates 
for further analysis. This yielded a final list of 108 iden- 
tified putative DML that are associated with 63 different 
DMRs from 62 genes (one gene with both paternal and 
maternal DML) (Table S2 in Additional file 2). Of the 
63 DMRs, 37 are maternally methylated and 26 are 
paternally methylated (Figure 3). These imprinted DMRs 
are distributed across the whole genome, with chromo- 
some 7 containing the highest number (nine DMRs), 
while chromosomes 13, 21 and Y are the only chromo- 
somes for which no DMRs were identified (Figure 3). 

As copy number variation (CNV) can be a potential 
bias for methylation [26], we referred to the UCSC Gen- 
ome Brower database (hgl8) (http://www.genome.ucsc. 
edu) and found that the locations of 37 of the 108 
probes overlap with known CNVs (Table S2 in Addi- 
tional file 2). However, any effect of the CNVs on 
methylation of the candidate sites identified by our cri- 
teria was minimal, since the methylation of maternal 
and paternal DML were clearly separated from each 
other without much overlap (Figure SID in Additional 
file 1). Similarly, differences between the two groups are 
unlikely to be caused by differences in genetic sequence 
polymorphisms that influence methylation, as this would 
require all ten diandric placentas, by chance, to be of a 
different genotype from all ten dygynic placentas. 

Validation of DNA methylation patterns of identified 
putative imprinted DMRs 

The microarray included 374 CpG sites in the promoter 
regions of 59 genes that have previously been reported 
to be imprinted in humans based on the literature [12] 
and information in Internet databases (http://igc.otago. 
ac.nz/ and http://www.geneimprint.com/) (Table S3 in 
Additional file 2). For nine of these genes {PRIM2A, 
IGF2R, TFP12, C0PG2, KLF14, ABCAl, INPP5F, 
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Figure 2 Scatterplots of average methylation of paternal and maternal DML. (A) and (C) Average methylation values in normal placentas 
(x-axis) plotted against digynic triploids (Avg G), diandric triploids (Avg A) and CHMs (Avg C) show high correlation. (B) and (D) Average 
methylation values in CHMs (x-axis) plotted against digynic triploids (Avg G), diandric triploids (Avg A) and normal placentas (Avg N). 



IGF2AS and BLCAP), the included CpG sites were 
unmethylated (or very lowly methylated) in the normal 
human placentas, as well as in the triploid placentas 
(Table S3 in Additional file 2). The majority (39 of 50) 
of the remaining genes showed differences between dia- 
ndric and digynic triploids at one or more of the asso- 
ciated CpG sites (^-test, P < 0.05), though typically not 
all CpG sites were differentially methylated. 

Among the 62 genes identified with parent-of-origin- 
dependent DMRs (using the stricter criteria of <0.1% 
FDR and absolute average mean difference >15%), 18 
are known imprinted genes associated with 15 distinct 
DMRs. Two of the identified DMRs, associated with the 
imprinted genes CDKNIC and RASGRFl, have been 
reported only in mice and not in humans [3,12] (Table 
1). While our strict selection criteria yielded only 18 of 
the 39 known imprinted genes that were statistically sig- 
nificantly different between diandric and digynic triploi- 
dies using an uncorrected P < 0.05 (Table S3 in 
Additional file 2), the missed cases were largely due to 



the mean difference being less than the 15% average 
methylation difference cutoff. Eleven of the fifteen 
imprinted DMRs are known to be ICRs with a parental 
origin of methylation concordant with what we observed 
based on the comparison of triploidies (Table 1). 

To confirm the methylation differences using an inde- 
pendent approach, we performed bisulphite pyrosequen- 
cing for a subset of the novel imprinted DMRs. For this 
purpose, ten DMRs were selected on the basis of their 
low FDR {FAM50B, MCCCl, DNAJC6, SORD and 
RH0BTB3) or their biological significance to the pla- 
centa {APC, DNMTl, IGFBPl, LEP and RASGRFl), A 
high correlation between the values obtained by micro- 
array and pyrosequencing was observed (r = 0.85 to 
0.98; P < 0.0001) (Figures S2A to S2J in Additional file 
1). Specifically, the DNA methylation patterns observed 
by pyrosequencing were concordant with those found 
by microarray for both (1) CpG sites analyzed by micro- 
array and their the proximal CpG sites within the pyro- 
sequencing assays (Figures S3A to S3J in Additional file 



Yuen et al. Epigenetics & Chromatin 201 1, 4:10 
http://www.epigeneticsandchromatin.eom/content/4/1/10 



Page 6 of 16 



Maternal DMR: 37 
Paternal DMR: 26 





<— C3orf62 


















S^NAPILS 






















^ MCCCl 
'*~FGF12 







. FAM50B 
' CD83 



RH0BTB3 
I NUDT12 

~ ARC FU 37396 



IGFBPl 
FIGNLl 



LEP 
PARP12 ■ 



H! 



' SGCE 
MEST 



hi 



_IRF7 
~ KCNQ1 



-SORD 
- RASGRFl 



ZNF597 ■ 

TBX6 ■ 



AKAPIO 

ccRio — ►r 



REEP6 — ►r— ^ ^ 



GNAS 

SAMDIO ' 



. L3MBTL 
- GNAS 




Figure 3 Location of the 63 identified DMRs in the genome. Relative location of the identified 37 maternal DMRs and 26 paternal DMRs are 
shown in the human genome according to the genomic sequence released in 2006 in the UCSC Genome Browser database (hgl8). Paternal 
DMRs are highlighted in blue and maternal DMRs are highlighted in red. Known imprinted genes are boldfaced and underlined. Chromosome 7 
contains the highest number of DMRs (nine DMRs), while there are no DMRs identified on chromosomes 13, 21 and Y. 



1) and (2) the average methylation levels of all CpG sites 
covered by pyrosequencing (Figures S4A to S4J in Addi- 
tional file 1). DNA methylation levels of the selected 
loci were also assessed in sperm DNA and all were 
unmethylated (data not shown), suggesting they may be 
either secondary DMRs or maternal imprinted DMRs. 

We further evaluated DNA methylation for two genes, 
FAM50B and MCCCl, which contain SNPs with high 
average heterozygosity (about 0.4) in the proximal pro- 
moter regions that can be used to distinguish alleles 
(Figures 4 A and 4F). Most of the other identified genes 
do not contain common SNPs in the nearby analyzed 
regions that could be used for this purpose. Bisulphite 
cloning and sequencing confirmed monoallelic methyla- 
tion patterns for both DMRs (Figures 4C and 4H) and 
maternal origin of allelic methylation that was concor- 
dant with that predicted by the triploidy comparison 
(Figures 4B and 4G). Furthermore, allelic expression 
analysis showed preferential expression of the unmethy- 
lated paternal allele at the proximal promoter regions 
(Figures 4E and 41), which is consistent with an inverse 
correlation between methylation and expression. As 



allelic methylation can occur in a SNP-dependent man- 
ner [27], we developed a methylation-specific pyrose- 
quencing assay for FAM50B to evaluate allelic 
methylation in additional samples. This same approach 
could not be applied to MCCCl, because its interro- 
gated SNP is located at a CpG site. The results of the 
FAM50B assay were concordant with cloning and 
sequencing results for the same placental sample (Fig- 
ures 4C and 4D). As methylation was found in associa- 
tion with either allele (A or G at rs2239713) among 12 
heterozygous normal term placental samples and ten 
heterozygous maternal blood samples (Table S4 in Addi- 
tional file 2), the allelic methylation is not linked to the 
SNP genotypes, at least for this DMR. 

Since diandric triploid placentas tend to be associated 
with trophoblast hyperplasia [17], we considered the 
possibility that some of the identified imprinted DMRs 
reflect differences in methylation between the tropho- 
blast and mesenchyme, the two components of the 
chorionic villi [28]. To address this hypothesis, we used 
a nonimprinted, trophoblast-specific unmethylated 
region, EDNRB (Figure S5A in Additional file 1) to 
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^DMR, differentially methylated region; ICR, imprinting control region; "^tissue- 
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compare the methylation levels between diandric and 
digynic triploid placentas. However, we did not find a 
difference in methylation levels between them at this 
site (Figure S5B in Additional file 1). Likewise, we did 
not find any differences in allelic methylation between 
trophoblast and mesenchyme for the novel identified 
imprinted gene MCCCl (Figures S5C and S5D in Addi- 
tional file 1). However, DNAJC6 and RASGRFl showed 
differential methylation between trophoblast and 
mesenchyme (Figure S6 in Additional file 1), which may 
represent cell-type-specific imprinting. 

Confirmation of parent-of-origin allelic expression for the 
identified putative imprinted genes 

While the existence of an imprinted DMR is thought to 
be predictive of imprinting at the gene expression level, 
proving this is complicated by the fact that imprinted 
DMRs may exist in association with imprinted genes 
even in tissues where the gene is not expressed or is 
expressed in a biallelic manner. Furthermore, to be 
informative for demonstrating monoallelic expression, 
the placental sample must be heterozygous for an 
expressed SNP. To be informative for parental origin of 
the expressed allele, the mother must additionally be 
homozygous for the same SNP. 

As many previously reported imprinted genes are 
expressed in an imprinted manner in placenta, and 
since we needed to screen many placenta-mother pairs 



to find informative cases, we proceeded to investigate 
the parental origin of allelic expression for the novel 
putative imprinted genes using a high-throughput geno- 
typing approach, specifically the iPLEX Gold assay on 
the MassARRAY platform (Sequenom, Inc., San Diego, 
CA, USA). We selected 38 of 45 genes associated with 
novel imprinted DMRs (the 45 putative imprinted genes 
including RASGRFl, for which imprinted expression has 
not been reported in humans) on the basis of the avail- 
ability of an exonic SNP with high average heterozygos- 
ity (>0.1) and the presence of expression in the placenta 
according to the GNF atlas database (http://biogps.gnf. 
org/). In addition, two exonic SNPs from 1GF2 were 
included as positive controls, since 1GF2 is well known 
to exhibit imprinted expression in human placentas. 
Thus, a total of 40 SNPs were genotyped in 27 mater- 
nal-foetal pairs, including DNA from maternal blood 
and the corresponding foetal normal term placenta, as 
well as cDNA from the same placenta. 

Of these 40 SNPs, 7 did not pass the quality control 
criteria (<70% calls or presence of severe allelic bias) 
and 3 had no informative (heterozygous) genotypes in 
foetal DNA, leaving a total of 30 SNPs for analysis 
(Table S5 in Additional file 2). The two SNPs from 
IGF2 showed the expected paternal allelic expression in 
all informative cases (Table S5 in Additional file 2). Of 
the 28 novel putative imprinted genes, 11 showed 
monoallelic expression in at least a portion of informa- 
tive samples (Table 2). Among these 11 genes, 8 had 
cases informative (homozygous) in maternal blood for 
parental origin assessment. Since most CpG sites in the 
microarray are located at the proximal promoter regions 
of the genes, we assumed that the DNA methylation 
would most likely correlate with silencing for all these 
genes. Six genes {FAM50B, DNMTl, RH0BTB3, 
ARMC3, AIFM2 and LEP) showed parent-of-origin- 
dependent expression that matched that predicted by 
the parental origin of the DMRs, while two others 
{MOVIOLI and ST8SIA1) showed parental expression 
opposite that predicted in one or more informative 
cases (Table 2). For FAMSOB and RH0BTB3, monoalle- 
lic expression for both reciprocal forms of the SNP was 
also observed. Some genes with imprinted DMRs may 
not show allele-specific expression biases because of the 
presence of tissue-specific or gestational age-specific 
imprinting that is further regulated by DNA methylation 
at other nearby sites. 

A number of genes did not consistently show monoal- 
lelic expression using the iPLEX Gold assay. For exam- 
ple, for LEP, only 1 of 15 samples was scored as 
monoallelic using this approach. To evaluate the sensi- 
tivity of the iPLEX Gold genotyping assay for detecting 
allelic biases in expression, we developed an RNA-speci- 
fic genotyping pyrosequencing assay for LEP, Although 
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Figure 4 Identification of imprinted DMRs at the proximal promoter regions of FAM50B and MCCCl. (A) and (F) Schematics showing the 
positions of methylation assays (Biseq: bisulphite cloning and sequencing assay; eg code: probe number of lllumina assay; and Pyro: bisulphite 
pyrosequencing assay) and SNP locations relative to the genes. Arrow directions represent the transcriptional directions for the genes. Genomic 
coordinates were retrieved from the UCSC Genome Brower database (hgl8). (B) and (G) Box plots showing the methylation levels of samples 
from each placental group for the DMRs analyzed by bisulphite pyrosequencing. Both DMRs in FAM50B and MCCCl have higher methylation in 
digynic than diandric triploid placentas, while they have intermediate methylation in normal placentas and particularly low methylation in CHMs. 
(C) and (H) Bisulphite cloning and sequencing showing parental origins of methylated and unmethylated alleles (M: maternal alleles; P: paternal 
alleles). Parental origin was determined by genotyping heterozygous informative SNPs for each sample. The DMRs in both FAM50B and MCCCl 
are maternally methylated. Each black circle represents a methylated CpG dinucleotide, and each white circle represents an unmethylated CpG 
dinucleotide. (D) Quantitative genotyping of methylated alleles by pyrosequencing. SNP rs2239713 is homozygous (GG) in maternal DNA and 
heterozygous (AG) in foetal (placental) DNA (dispensation order: AAG). Genotyping of the placental sample using a methylation-specific 
pyrosequencing primer shows a homozygous (GG) pattern, indicating that the DMR associated with the maternally inherited 'G' allele is 
methylated while the one associated with the paternal 'A' allele is not. (E) and (I) Quantitative genotyping of expressed alleles by 
pyrosequencing. Both SNPs (E) rs6597007 (dispensation order: GGC) and (I) rs937652 (dispensation order for DNA genotyping: CG; dispensation 
order for RNA genotyping: CCG) are homozygous in maternal DNA and heterozygous in foetal DNA. Genotyping of cDNA shows a bias towards 
preferential expression of the paternal alleles. *The pyrosequencing primers used for cDNA genotyping (intron-spanning) in MCCCl were 
different from those used for DNA genotyping (Table S10 in Additional file 2), so the peak ratio shown in genotyping the pyrogram of cDNA 
does not correspond to that for DNA. 



the two methods were correlated (r = 0.64; P < 0.02), we 
found that pyrosequencing was more likely to detect 
preferential allelic expression, with 5 of 12 informative 
cases exhibiting a <0.3 allelic ratio by pyrosequencing 
(Table S6 in Additional file 2). Furthermore, in case 
PM155 for MCCCl, we found preferential paternal 



allelic expression by pyrosequencing (Figure 41), but not 
by iPLEX Gold genotyping (Table 2). Thus, the iPLEX 
Gold assay may not be sufficiently sensitive to detect 
more subtle allelic expression bias, that is, in circum- 
stances where there is a mix of cells with biallelic and 
monoallelic expression. 
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Table 2 Eleven genes associated with candidate imprinted OMRs with confirmed monoallelic expression^ 

Gene DMR SNP Monoallelic expression, Monoallelic expression observed for Matched expected parental origin, 

observed/total (%) reciprocal SNP^ observed/total 



FAM50B 


M 


rs6597007 


9/9 (100) 


Y 


5/5 (100) 


DNMTl 


M 


rsl 6999593 


1/1 (100) 




1/1 (100) 


MOVIOLI 


P 


rs96 17066 


8/9 (89) 


N 


1/3 (33) 


RH0BTB3 


M 


rs34896 


3/4 (75) 


Y 


2/2 (100) 


SNCB 


M 


rs2075667 


3/4 (75) 


N 


Nl 


ARMC3 


M 


rsl 2259839 


2/3 (67) 


N 


2/2 (100) 


ST8SIA1 


M 


rs4762737 


2/3 (67) 


Y 


0/1 (0) 


ARHGAP4 


P 


rs2070097 


1/2 (50) 




Nl 


AIFM2 


M 


rs7908957 


2/8 (25) 


N 


1/1 (100) 


MCCCl 


M 


rs937652 


2/8 (25) 


Y 


Nl 


LEP 


P 


rs2 167270 


1/15 (7) 




1/1 (100) 



^DMR, differentially methylated region; SNP, single-nucleotide polymorphism; Nl, not informative. "^Where both alleles of SNP were observed to be expressed 
among cases with monoallelic expression; this is impossible to determine if only one case showed monoallelic expression. ^Number of cases matching the 
expected parental origin of those cases informative with regard to determining parent of origin. 



Tissue-specific and gestational age-specific methylation of 
imprinted OMRs 

To study tissue-specific effects and the effect of gesta- 
tional age on methylation of the putative imprinted 
DMRs, we further compared methylation at these sites 
among three types of foetal somatic tissues (eight brain 
samples, twelve kidney samples and eleven muscle sam- 
ples) and two sets of placentas with different gestational 
ages (ten midgestation and ten term placentas) that had 
been run in the same Infinium HumanMethylation27 
BeadChip array. 

For tissue-specific methylation analysis, we compared 
the DNA methylation levels of the 108 DML (probes) 
associated with the 63 imprinted DMRs in five tissues 
(brain, kidney, muscle, midgestation placenta and blood). 
Multiclass comparison from SAM was performed with 
1,000 permutations. Using a <0.1% FDR cutoff, 53 probes 
of 46 imprinted DMRs showed differential DNA methyla- 
tion between tissues (Table 3 and Table S7 in Additional 
file 2). Placenta-specific methylation was observed for 31 
of these probes (26 imprinted DMRs), with the average 
methylation being more than 15% higher in placenta than 
in any other tissues (Table 3 and Table S7 in Additional 
file 2). A change in methylation of placenta by gestational 
age was found for 12 probes from ten DMRs using the 
same statistical criterion (<0.1% FDR) (Table 3 and Table 
S8 in Additional file 2). Thus, imprinted DMRs can show 
both tissue-specific and gestational age-specific DNA 
methylation. Nonetheless, 14 of the imprinted DMRs 
showed constant methylation between different tissues 
and gestational ages (Table 3 and Table S9 in Additional 
file 2), 11 of which are in ICRs from known imprinted 
genes. Three identified imprinted DMRs associated with 
FAM50B, FGF12 and IRF7 also remained constant across 
samples and are thus potential ICRs or primary DMRs. 



The complexity of DNA methylation associated with 
imprinted genes can be illustrated by the data for three 
known imprinted genes, GNAS, CDKNIC and MEST, 
for which multiple probes were present on the Infinium 
HumanMethylation27 BeadChip array. For GNAS, the 
array contains probes for 30 CpG sites mapping across 
three promoter regions of three alternative transcripts 
{NESPSS, GNASXL and exon lA of GNAS) (Figure 5A). 
As previously reported, the paternal DMR is located at 
the promoter of NESPS5 transcript (Figure 5B), while 
the maternal DMR is located at the promoter of 
GNASXL [29]. While most of the CpG sites have more 
or less equal average methylation across the locus, 
cgl5160445 to cgl683351 and cg01565918 show clear 
tissue-specific methylation across different tissues (Fig- 
ures 5B to 5D). For CDKNIC, there are eight probes 
present in the array (Figure 5E). A previously unidenti- 
fied paternal DMR was identified at the promoter region 
of this gene through our comparison of triploids (Figure 
5F). Interestingly, not only is the imprinted DMR itself 
tissue-specific (that is, it is a secondary DMR) (Table 3) 
but there is also a probe (cg209 19799) that shows differ- 
ential methylation across different gestational ages 
(Figure 5G) and tissues (Figure 5H). Likewise, for 
MEST, for which ten probes span two regions of the 
gene (Figure S7A in Additional file 1), an imprinted 
DMR can be found in one region (Figures S7B and S7C 
in Additional file 1), while tissue-specific and gestational 
age-specific methylation is observed in another region of 
the MEST promoter (Figures S7C to S7G in Additional 
file 1). 

Discussion 

Many efforts have been made to identify imprinted genes 
in the human genome because of their importance in 
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Table 3 DNA methylation of identified DMRs in different tissues and gestational ages^ 



Index 


Gene 


Chromosome 


Tissue-specif ic'^ 


Change in gestation^ 


Stable non-tissue-specific*^ 


Known imprir 


] 


DNA IC6 


] 




Y 


N 


N 


2 


LASS2 


1 


Y^ 


Y 


N 


N 


3 


PEX5 


1 2 


Y^ 


Y 


N 


N 


4 


RASGRFl 


1 5 


/ 


N 


N 


N 


5 


AKAPin 


1 7 


Y^ 


N 


N 


N 


6 


AIFM2 


10 




N 


N 


N 


7 


A PC 


5 


/ 


N 


N 


N 


3 


ARHGAP4 


X 


/ 


N 


N 


N 


9 


ARMC3 


10 


Y^ 


N 


N 


N 


10 


C3orf62 


3 


Y^ 


N 


N 


N 


1 1 


CD83 


6 


Y^ 


N 


N 


N 


12 


CMTM3 


16 


/ 


N 


N 


N 


1 3 


DNMTl 


1 9 


Y^ 


N 


N 


N 


14 


G0S2 






N 


N 


N 


1 5 


GATA4 


8 


/ 


N 


N 


N 


1 6 


LEP 


7 


/ 


N 


N 


N 


1 7 


MCCCl 


3 


Y^ 


N 


N 


N 


1 8 


NlJDTl? 


5 


Y^ 


N 


N 


N 


1 9 


PCK2 


14 


Y^ 


N 


N 


N 


20 




5 


/ 


N 


N 


N 


21 


SLC46A2 


9 


Y^ 


N 


N 


N 


22 


SNCB 


5 




N 


N 


N 


23 


SORD 


1 5 


/ 


N 


N 


N 


24 


ST85IA 1 


1 2 


/ 


N 


N 


N 


25 
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Y^ 


N 


N 


N 


26 


TMFM 1 7 

1 IVILIVI 1 / 


2 


Y^ 


N 


N 


N 


27 
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1 7 


Y^ 


N 


N 
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N 


N 


N 


29 




6 


Y 
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N 
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30 


D IRAS 3 
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32 
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Y 
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33 




1 1 


Y 


N 


N 


Y 


34 
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1 1 


Y 


N 


N 


Y 


35 
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36 
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37 


PEGIO 


7 


Y 


N 


N 


Y 


38 


GlOnrfVS 


10 


Y 


N 


N 


N 


39 


CCRIO 


1 7 


Y 


N 


N 


N 


40 


CYP2W1 


7 


Y 


N 


N 


N 


41 


FIGNL 1 


7 


Y 


N 


N 


N 


42 


IGFBPl 


7 


Y 


N 


N 


N 




hAr)\/WI 7 
IVIkJV I UL 1 


zz 


Y 
Y 


M 
In 


M 
In 


M 
In 


44 


P2RY6 


11 


Y 


N 


N 


N 


45 


PARP12 


7 


Y 


N 


N 


N 


46 
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47 
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48 
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Table 3 DNA methylation of identified DMRs in different tissues and gestational ages^ (Continued) 



52 


(oHd lU 


7 


M 

IN 


M 

IN 


V 

Y 


V 

Y 


JO 


ML J 1 


/ 


M 


M 

In 


V 

Y 


V 

Y 


54 


NAriLD 


4 


N 


N 


Y 


Y 


55 


PEGS 


19 


N 


N 


Y 


Y 


56 


PLAGLl 


6 


N 


N 


Y 


Y 


57 


SGCE 


7 


N 


N 


Y 


Y 


58 


SNURF 


15 


N 


N 


Y 


Y 


59 


ZIM2 


19 


N 


N 


Y 


Y 


60 


ZNF597 


16 


N 


N 


Y 


Y 


61 


FAM50B 


6 


N 


N 


Y 


N 


62 


FGF12 


3 


N 


N 


Y 


N 


63 


IRF7 


11 


N 


N 


Y 


N 



^DMR, differentially methylated region; "^multiclass comparison of methylation level in brain, kidney, muscle, midgestation placenta and blood with FDR <0.1%; 
'^Multiclass comparison of methylation level in early-gestation, midgestation and term placenta, FDR <0.1%; ^DMRs with no statistically significant changes in 
methylation level in different tissues and gestational ages; ^Based on the public databases (http://igc.otago.ac.nz/ and http://www.geneimprint.com/); ^placenta- 
specific methylation. 



foetal growth and development and their potential for 
dysregulation [11,12]. Most known imprinted genes to 
date were first identified in mice, but many imprinted 
genes are not conserved across species [5]. In the present 
study, we utilized diandric and digynic triploid placentas 
to map imprinted DMRs, sites that are typically asso- 
ciated with imprinted genes, in the human genome. We 
identified 11 of the 18 previously reported human ICRs 
covered by the lUumina Infinium HumanMethylation27 
BeadChip panel, with additional ones which showed dif- 
ferences that were insufficient to reach our stringent sta- 
tistical criteria. Furthermore, we confirmed the parent-of- 
origin dependence of methylation and expression in a 
subset of our candidate novel imprinted genes on the 
basis of independent experiments. 

This approach improves upon previous strategies for 
mapping imprinted DMRs, such as comparing partheno- 
genotes (ovarian teratomas) and androgenotes (CHMs) 
[13,14], which is limited by the grossly abnormal nature 
of these samples, or comparing maternal and paternal 
uniparental disomies (UPDs) [30,31], which is restricted 
by the rarity of UPDs for many chromosomes and the 
limited tissues available for analysis. Although triploid 
placentas do exhibit some abnormal pathology, their cel- 
lular composition is comparable and methylation pro- 
files of both types of triploidy were closely correlated 
with chromosomally normal placentas (r = 0.99). In 
comparison, a previous study showed that mature ovar- 
ian teratomas have a methylation profile more similar to 
that of blood (r = 0.94) than to either CHMs (r = 0.84) 
or normal placentas (r = 0.88) [14]. Genome-wide tran- 
scriptome analysis has also been used to identify 
imprinted genes [11,32], but it is gene expression- and 
SNP-dependent; thus, imprinted genes with tissue-speci- 
fic expression or lacking a heterozygous exonic SNP 
would be missed. 



As demonstrated, tissue-specific methylation of 
imprinted DMRs or their flanking regions can readily be 
assessed by comparing methylation profiles of a variety 
of tissues, allowing a comprehensive analysis of tissue- 
specific methylation regulation at complex loci, such as 
GNAS [29]. The regional dependent methylation pat- 
terns in the promoters of imprinted genes show the 
importance of locating the specific CpG sites defining 
the imprinted DMRs when studying the dynamics of 
promoter DNA methylation at such genes. While in the 
present study we identified only loci that demonstrated 
parent-of-origin-dependent differential DNA methyla- 
tion in placenta, most known imprinted genes show par- 
ent-of-origin-specific expression in this organ [2]. 
Furthermore, as diandric and digynic triploids can both 
exist as foetuses, additional comparisons could be used 
to identify any potential genes that exhibit imprinting 
specifically in other tissues. This study is limited by the 
low coverage of CpG sites in the array (about two CpG 
sites on average for each proximal promoter region of 
genes), which reduces its power to identify imprinted 
DMRs as these may be limited to specific regions within 
the promoter. This analysis could thus be extended 
further by using microarray or whole-genome sequen- 
cing with greater coverage of the genome. 

Overall, the number of imprinted DMRs identified in 
the present study was less than that predicted by bioin- 
formatics approaches [33]. However, the stringent selec- 
tion criteria (<0.1% FDR and absolute average 
methylation difference >15%) that we used to pick the 
top candidate sites caused an underestimation of the 
number of imprinted loci. Many more candidate 
imprinted DMRs can be identified with this data set 
using lower thresholds. In fact, many known imprinted 
genes that we failed to identify on the basis of these cri- 
teria did show nominally significant {P < 0.05 without 



Yuen et al. Epigenetics & Chromatin 201 1, 4:10 
http://www.epigeneticsandchromatin.eom/content/4/1/10 



Page 12 of 16 



Scale: 
Chr20: 



GNAS B 
Region 1 > 



NESP55 
GNASXL 
Exon 1A 

CpG 
Islands: 



56855000 I 



GNAS Region 2 I 



20 kb I 

56865000 I 



56875000 I 



56885000 I 



GNAS Region 3 



56895000 I 



GNAS Region 1 ^ 



GNAS Region 2 



GNAS Region 3 



t 0.4 



0.8 



0.6 - 



1 n 



0.6 



0.2 




0.4 



— I — I — I — I — I — I — I — I — I — I — -1 — -1— 




-Diandric 
-Digynic 
- Normal 
-CHM 



-PLN(M) 

-MUS 

-BRN 

-KID 

-WB 



LO 00 O CM CO O) 

o o) ^ a> 00 00 

00 05 CO 00 00 h~ a> 

CO 1^ 00 Lo a> LO CNj 

OJ 0> LO 1^ T- LO 00 

uo Lo Lo a> CD 00 LO 



■ 00 ^ 1^ CD CD 



00 CD 
00 00 CD 
CD t- - 



CM °0 ? 

C^CDSSLOIZOOh- 

^^SR;$;<M?5°ooo 

05co°ooocdOc\I'^ 



0)0)0)0 



ocNiT-ocNJi-frrrCNj 

0)0)0)0)0)9i'S^C" 



00 00 00 o T- Lo o 

T- CNJ O 00 CN CM O 

CD CD 00 CO 00 CM C35 

uo T- CO uo 00 -"^l- 

CD T- CNI 00 00 CD 

UO O O CD 00 LO o 

T- O LO CM CD CD 

O T- CN CM CN O O 
0)0)0)0)0)0)0) 



Scale: 
Chr11: 



CDKN1C 
CpG 



1 kb I- 



cg031 43849 



2862000 I 
cg23146358 



2863000 I 



2864000 I 



cg04548378 



Icg20919799 I 

cg23548479 



cg1 77381 33 I 
cg25215340 I 



cg05559445 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0 




-Diandric 
-Digynic 
Normal 
-CHM 



0.7 
0.6 - 
0.5 
0.4 
0.3 
0.2 
0.1 
0 



-PLN(E) 
-PLN(M) 
-PLN(T) 




-PLN(M) 

-MUS 

-BRN 

-KID 

-WB 



CDOOOOCDCDOOO* 

'itior^cDi^oo'itS 
ooooooi^'^T-oo^ 

COCDOOCDOOOOUOSC 
'^^'^■<-^00t-O5 

T-T-LOCDLOI^CN^ 

ooco-^ocoi^LO"^ 

OC\IOC\ICNjT-CNj"5 
0)0)0)0)0)0)0) 
OOOOOOOC3) 



CDOOOOCDCDOOO* 

'itior^cDr^oo';!-^ 
ooooooi^'^T-oo^ 
oocdoocdooooloSl 

^^^t-tI-oOt-CD 

T-T-LocDuni^c\i"5 
00 00-^ocoi^un"^ 

OC\IOC\IC\l-<-C\l"5 
0)0)0)0)0)0)0) Si 

ooooooqc^ 



CDOOOOCDCDOOO* 

■^Lor^cDr^oo-^iC 
ooooooi^'itT-003 
cocdoocdoooolo^ 

ooooti-oooi^lo"^ 

OC\IOC\IC\It-CNJ^ 
0)0)0)0)0)0)0) 

ooooooog> 



Figure 5 Illustration of tissue-specific and gestational age-specific methylation at the proximal promoter regions of GNAS and 

CDKN1C. (A) and (E) Schematics showing the positions of the lllumina Infinium probes relative to the genes and transcripts. Arrow directions 
represent the transcriptional directions for the genes or transcripts. Genomic coordinates were retrieved from the UCSC Genome Brower 
database (hgl8). (B) through (D) Average methylation levels of the lllumina Infinium probes in different placental groups (top) and in different 
tissues (bottom). Probe numbers are shown on the x-axes in the bottom panels divided into (B) GNAS region 1, (C) GNAS region 2 and (D) GNAS 
region 3 according to their proximity to the known transcripts. Tissue-specific methylation can be found from cgl 51 60445 to eg 16833551 in 
GNAS region 2 and at cg01565918 in GNAS region 3. (F) through (H) Average methylation levels of the lllumina Infinium probes of CDKNIC in 
(F) different placental groups, (G) different gestational ages of placenta and (H) different tissues. Probe numbers are shown on the x-axes. Both 
tissue-specific and gestational age-specific methylation can be found at cg20919799. PLN(E): early gestation placenta; PLN(M): midgestation 
placenta; PLN(T): term placenta; MUS: muscle; BRN: brain; KID: kidney; WB: whole blood. 
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correction for multiple comparisons) DNA methylation 
differences between diandric and digynic triploids (Table 
S3 in Additional file 2). For instance, a recently con- 
firmed imprinted gene, RBI [34], was significantly differ- 
entially methylated between diandric and digynic 
triploidies (<0.1% FDR), with a methylation pattern con- 
sistent with that of a maternal DMR (data not shown). 
However, it was excluded because its absolute average 
methylation difference between diandries and digynies 
was only 14%. While we expect imprinted DMRs to 
show a difference of 33% between the triploid groups, 
smaller differences may be observed owing to a lack of 
complete methylation at all CpG sites on the inactive 
allele or to the presence of a mix of cell types, only 
some of which are imprinted. Similarly, methylation on 
the inactive X chromosome in females is incomplete 
(much less than 50%) in the placenta for gene promoter 
regions that are typically methylated at 50% in somatic 
tissues, despite still showing a significant increase in 
methylation relative to male placenta [22]. 

Only some of the novel putative imprinted DMRs 
could be confirmed to show monoallelic expression, and 
others did not show strict parent-of-origin expression in 
all placentas (Table S5 in Additional file 2). There are 
several possible explanations. First, there may be cell- or 
tissue-specific imprinting confounding the ability to 
detect a difference in whole villous samples from term 
placentas. Many known imprinted genes show imprinted 
expression only in specific placental cell types, for exam- 
ple, Mash2 in mice, which is differentially expressed 
only in diploid trophoblast cells of the postimplantation 
embryo [35], and STOXl in humans, which is mater- 
nally expressed in extravillous trophoblast cells [36]. 
Given the highly heterogeneous cell types present in the 
placenta [28], nonimprinted expression in some cells 
may mask allelic expression in others. The possibility 
that cell heterogeneity exists for the DMRs identified in 
the present study is supported by the observations that 
(1) average methylation of some candidate DMRs was 
not the expected 50% in normal placentas (Figure S3 
in Additional file 1) and (2) DNAJC6 and RASGRFl 
showed differential methylation between trophoblast 
and mesenchyme (Figure S6 in Additional file 1). Sec- 
ond, as we have shown, the iPLEX Gold assay may not 
be sensitive enough to pick up subtle allelic expression 
biases (Table S6 in Additional file 2). 

Third, there may be alternative transcripts regulated 
by alternative promoters that are not imprinted, so the 
observed expressed allelic ratio at particular SNPs may 
be complicated by the synergic effect of multiple tran- 
scripts. Such complex regulation is observed for known 
imprinted genes such as GNAS, CDKNIC and MEST 
(Figure 5 and Figure S5 in Additional file 1). However, 
allelic expression from either parent in some genes. 



such as MOVIOLI and ST8SIA1, suggests that some of 
the identified DMRs may be random monoallelically 
expressed genes instead of imprinted genes (with DNA 
methylation differences between diandries and digynies 
occurring by chance). 

The validation of all the putative imprinted DMRs we 
identified is limited by the number of samples and com- 
mon SNPs within regions and by the availability of 
intact mRNA from the pathological specimens. A proper 
validation experiment to demonstrate that the DMRs we 
have identified are associated with imprinted methyla- 
tion and gene expression requires being able to trace 
the parental origin of the methylated and expressed 
alleles in multiple members of the same family, which 
can be done in mice but is impractical and ethically 
impossible to do across multiple tissues in humans [37]. 
The best alternative is to trace the origin of the methy- 
lated allele and the expressed allele in multiple indivi- 
duals. This requires a SNP adjacent to the methylation 
site that is heterozygous in the test sample but homozy- 
gous in one parent. Using this strategy, we demon- 
strated for FAM50B that (1) a maternal origin of the 
methylated allele in placenta and blood from multiple 
individuals and on reciprocal genetic backgrounds and 
(2) the paternal allele is expressed with either SNP allele 
in the placenta, thus ruling out the possibility of a 
genetic effect. Confirming that an imprint represents a 
primary imprinted DMR requires detailed investigations 
of postfertilization imprinting dynamics, which is diffi- 
cult to perform in humans. Nonetheless, we showed 
that the methylation level of FAMSOB is similar in mul- 
tiple tissues and is unmethylated in sperm, suggesting 
that it is likely to be a primary maternal DMR. During 
the revision of this manuscript, the maternal imprint of 
FAMSOB was also confirmed by other groups using 
similar validation methods [38,39]. The goal of the pre- 
sent study was to demonstrate the ability of our 
approach to identify imprinted DMRs, not to map and 
confirm every imprinted DMR on the array. Thus, the 
putative imprinted DMRs listed in the present study 
should be considered with caution, and further valida- 
tion is required. 

Two genes identified as potentially being imprinted in 
the present study, APC and DNMTl, were excluded as 
being imprinted in previous studies [40,41], while APC 
was reported as being imprinted in another study [42]. 
Of interest, DNMTl is a DNA methyltransferase that is 
important for the maintenance and establishment of 
DMRs in imprinted genes [43], while APC is a negative 
regulator of the Wnt signalling pathway, which has been 
implicated in the survival, differentiation and invasion of 
human trophoblasts [40]. Although Dnmtl was found to 
be dispensable for growth of the extraembryonic 
lineages in mice [44], it is not methylated at the 
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orthologous region in mice [41]. Both the APC and 
DNMTl DMRs were reported to be specifically methy- 
lated in primate placentas [45], suggesting that the 
potential imprinting marks of these genes emerged fairly 
recently in evolution. This is also consistent with the 
hypothesis that maternal imprints are under selective 
pressure during early development for methylation- 
dependent control [46]. This could occur by selecting 
genes with developmental advantage by gain of imprint- 
ing from epipolymorphisms [47]. 

Conclusions 

In conclusion, we have demonstrated that comparison of 
diandric and digynic triploids is an effective method for 
mapping imprinted DMRs in the human genome. This 
approach can be extended to different tissues, gestational 
ages or species, thereby generating a comprehensive view 
of imprinting regulation and evolution. The ability to 
map novel imprinted DMRs in the human genome 
should improve our understanding of the causes of pla- 
cental dysfunction and birth defects. With the rapid 
advancement of molecular genetics technologies, a com- 
plete map of imprinted DMRs may ultimately be gener- 
ated by the use of whole-genome sequencing. However, 
the present approach is a convenient, currently available 
and cost-effective method of imprinted gene mapping. 

Methods 

Sample collection 

This study was approved by the ethics committees of 
the University of British Columbia and the Children's & 
Women's Health Centre of British Columbia. Early 
gestation placental samples (ten diandric triploids, ten 
digynic triploids, six CHMs and ten normal controls) 
were obtained from spontaneous abortions examined in 
the Children's & Women's Health Centre of British 
Columbia pathology laboratory. The parental origin of 
triploids was determined by using microsatellite poly- 
morphisms as previously described [17-19], and these 
studies also allowed us to exclude maternal contamina- 
tion in the placental samples. Midgestation placental 
samples {n = 10) and foetal tissues (11 muscle samples, 
12 kidney samples and 8 brain samples) were obtained 
from anonymous, chromosomally normal, second-trime- 
ster elective terminations for medical reasons. Term pla- 
cental samples and the corresponding maternal blood 
samples were collected from Children's & Women's 
Health Centre of British Columbia with the women's 
written informed consent. For all placental samples, 
fragments of about 1 cm^ were dissected from the foetal 
side and whole villi were used for investigation. All tis- 
sues were karyotyped for chromosomal abnormalities, 
and genomic DNA was extracted from each tissue sam- 
ple using standard techniques. Total RNA was extracted 



from term placentas using an RNeasy kit (Qiagen, 
Valencia, CA, USA) according to the manufacturer's 
instructions. 

Illumina DNA methylation array 

Genomic DNA was bisulphite-converted using the EZ 
DNA Methylation Kit (Zymo Research, Orange, CA, 
USA) according to the manufacturer's instructions. Bisul- 
phite treatment converted unmethylated cytosines to ura- 
cils while leaving methylated cytosines unchanged. After 
DNA purification, bisulphite-converted DNA samples 
were randomly arrayed and subjected to the Infinium 
HumanMethylation27 BeadChip panel array-based assay. 
The array assays methylation levels at 27,578 CpG sites 
in the human genome. The methylation level for each 
CpG site was measured by the intensity of fluorescent 
signals corresponding to the methylated allele (Cy5) and 
the unmethylated allele (Cy3). Cy5 and Cy3 fluorescence 
intensities were corrected independently for background 
signal and normalized using GenomeStudio software 
(Illumina, Inc.). Continuous P values that range from 0 
(unmethylated) to 1 (methylated) were used to identify 
the percentage of methylation, from 0% to 100%, for each 
CpG site. The p value was calculated based on the ratio 
of methylated/ (methylated + unmethylated) signal out- 
puts. The detection P value of each probe was generated 
by comparison with a series of negative controls 
embedded in the assay. Probes with detection P values 
>0.05 in any of the samples were eliminated from the 
study. The correlation coefficient for technical replicates 
was >0.98. The microarray data from this study have 
been submitted to the NCBI Gene Expression Omnibus 
(http://www.ncbi.nlm.nih.gov/geo) under accession num- 
ber GSE25966. 

DNA methylation analyses for targeted loci 

Methylation-unbiased PCR and sequencing primers were 
designed based on the probe sequences provided by Illu- 
mina (Table SIO in Additional file 2). All primers were 
designed in regions free of known SNPs. Pyrosequencing 
was performed using a PyroMark MD system (Biotage, 
Uppsala, Sweden). The quantitative levels of methylation 
for each CpG dinucleotide were evaluated using Pyro 
Q-CpG software (Biotage). For bisulphite cloning and 
sequencing, the PCR product from individual samples 
was generated by using non-biotinylated primers (Table 
SIO in Additional file 2) and subsequently TA-cloned 
into the pGEM-T Easy Vector System (Promega, Madi- 
son, WI, USA). Individual clones were picked and PCR- 
amplified with SP6 and T7 promoter primers. PCR pro- 
ducts were sequenced by using Sanger sequencing. The 
sequencing data were analyzed using BiQ Analyzer Soft- 
ware [48], and sequences with less than an 80% bisul- 
phite conversion rate were eliminated from analysis. 
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SNP genotyping 

Multiplex genotyping of genomic DNA and cDNA was 
performed by using the iPLEX Gold assay on the Mas- 
sARRAY platform (Sequenom) at the Genome Quebec 
Innovation Centre (Montreal, PQ, Canada). Primers for 
SNP genotyping were designed by using primer design 
software from Sequenom (Table Sll in Additional file 
2). The primer extended products were analyzed and 
the genotypes were determined by mass spectrometric 
detection using the MassARRAY Compact System 
(Sequenom). Technical replicates showed a correlation 
of r = 0.92. Samples or SNPs with <70% conversion 
rates (calls) were eliminated. Genotyping by pyrosequen- 
cing was performed on a PyroMark MD System, and the 
relative levels of alleles for SNPs were evaluated by 
using PSQ 96MA SNP software (Biotage). Genotyping 
of exonic SNPs was carried out with cDNA prepared 
using either (1) the Omniscript Reverse Transcriptase 
Kit (Qiagen) followed by the iPLEX Gold assay or pyro- 
sequencing or (2) the Qiagen OneStep RT-PCR Kit fol- 
lowed by pyrosequencing. Primers for pyrosequencing 
genotyping were designed by using primer design soft- 
ware from Biotage (Table Sll in Additional file 2). PCR 
without reverse transcriptase was performed on each 
sample to confirm that there was no genomic DNA 
contamination. 

Statistical analysis 

Unsupervised hierarchical clustering of samples was 
done using Illumina GenomeStudio software. Differen- 
tially methylated probes in the Illumina Infinium 
HumanMethylation27 BeadChip array from each com- 
parison were identified using the siggenes package from 
R software with a cutoff of <0.1% FDR. FDRs were gen- 
erated after comparison of 1,000 random permutations 
between samples. The Pearson linear correlation coeffi- 
cient was used to determine the similarity of DNA 
methylation profiles between samples. The Database for 
Annotation, Visualization and Integrated Discovery 
(DAVID) program was used for gene ontology analysis 
using the total number of genes presented in the array 
as a background for comparison [49,50]. 

Additional material 



Additional file 1: Figures SI to S6. 
Additional file 2: Tables SI to Sll. 
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